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This report contains the results and details of a design study and 
analysis performed by Advanced Technology and Research, Inc. (ATR) for the NASA 
Goddard Space Flight Center (GFSC) as part of the overall Inertial Energy 
Storage Hardware Definition Study. The study was confined to include only the 
composite ring rotor element. The GSFC technical officer for this task was Mr. 
G. Ernest Rodriguez. 

The objectives of this study were: 

1. Detailed design definition of a composite ring rotor meeting the 
stated performance requirements. 

2. Definative specifications for ring rotor fabrication and assembly. 

3. Identification of potential sources for fabrication and assembly. 

This design effort centered around the use of graphite-epoxy composite 

materials. The study is significant because it resulted in the following: 

1. High energy density "thick ring" composite rotor design compatible 
with the "Mechanical Capacitor" concept. 

2. Unique multiring prestressed design configuration. 

3. Successful application of a modified Tsai-Hill Quadratic Interaction 
Analysis technique. 

4. Demonstrated usefulness and versatility of previously developed 
computer codes (FLYANS and FLYSIZE) , which include the presence of an 
iron inner ring. 


I. DESIGN BACKGROUND 


This report covers the design definition of the composite fiber ring-rotor 
portion of an Inertial Energy Storage Hardware Definition Study managed by the 
NASA Goddard Space Flight Center (GSFC). Genesis for this effort stems from 
early work done at GSFC on the "Mechanical Capacitor" program. This early work 
showed the potential merits of a magnetically suspended ring rotor with 
integral ironless armature motor/generator construction. The three major 
elements of this design concept are: 

1. Composite material ring-rotor flywheel 

2. Magnetic suspension 

3. Permanent magnet ironless armature motor/generator 

The GSFC selected item 1. as the logical starting point for a complete 
design definition of the inertial energy storage system. 

Stated requirements for this effort were to provide design definition and 
a specification for a Inertial Energy Storage wheel comprised of composite 
materials to yield high energy density and be capable of interfacing with an 
appropriate magnetic suspension system and motor/ generator. 


Performance requirements: 


Energy Storage 
(at max. speed) 

1.6 KWh 

Weight 

(total system) 

60 kg(132 lbm) 

Volume 

Less than 0.05 m3 (1.77 ft3) 

Environment 


Temperature 

-10°C to +50°C (14°F to 122°F) 

(operating) 

. -4 

Pressure 

(operating) 

< 10 torr 


The basic requirement that establishes the size of the storage element is 
the 1.6 KWh capacity. Before determining specific parameters it was necessary 
to select a best estimate composite material. Previous work supported by the 
Department of Energy, Mechanical Energy Storage Technology (MEST) program 
showed that graphite fiber epoxy matrix material demonstrated best overall 
technical performance and ease of manufacturing when compared to other 
materials and fabrication technique*,. Thus, a graphite/epoxy fiber (Celion 
6000) and an epoxy matrix (EPON 826) were selected for this study. 

Two form factors (ID/OD ratio) of 0.5 and 0.4 were selected and the design 
sizing and analysis were conducted using these stated parameters. 
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II. ROTOR ANALYSIS 
1. Introduction 

In designing a magnetically suspended flywheel rotor it is necessary to 
consider a reasonable geometry upon which the design/analysis will be based. 

♦ 

Early work in this area by Kirk, Studer and Evans of GSFC [1, 2, 3]* has shown 
that a pierced disk of uniform thickness provides a desirable rotor geometry 
from both a performance and manufacturing point of view. 

Shown in Figure 0.1 is a cross sectional view of the original GSFC 
magnetically suspended flywheel design. The original design consists of 2 
rings with the outermost ring being made of a filamentary wound composite 
material and the inner ring being made of continuous iron bonded to the 
filamentary wound ring. The stator of this design fits in the !! hele" of the 2 
ring rotor and it carries the magnetic suspension and motor/generator 
electronics. 

RCA [4] and Kirk [5], under contract to GSFC, have done additional work on 
the original GSFC design. Kirk concluded that a multiring, interference 
assembled rotor, such as shown in Figure H.2, would provide for substantial 
improvements in the original GSFC design. The modified GSFC design, shown in 
Figure U2, differs from the original GSFC design in the following areas: 

1. The rotor is composed of a number of individual filamentary wound rings, 
rather than being one continuous ling. 

2. The inside diameter (ID) to outside diameter (OD) rotor ratio (ID/0D) 
is smaller than the original GSFC design. 


^Brackets denote references at the back of this section. 





Outside Radius 



Cross sectional viev of the aodified GSFC design 
The design has a smaller ID/OD ratio and the iro 
ring is segmented. 




3. The innermost ring is msde of iron and is segmented into discrete pic 
shaped "chunks' 1 * 

Each of. the above changes was made in the origiual GSFC design in order to 
improve the overall performance of the energy storage system. The reasoning 
behind the changes has been documented by Kirk and Huntington [6, 7* 8, 9) and 
a brief explanation follows: 

1. The rotor is made of a number of composite material rinps which are 
interference assembled. The reason behind this change is to favorably 
prestress the rotor so higher rotational speeds and energy densities 
can be obtained before a limiting performance constraint is 
encountered. 

2. The ID/OD ratio has been lowered. The reason behind this change is 
that the original GSFC design was of a "thin hoop" type And suffered 
excessive "gap" growth between the rotor and the stator as it spun. 
Since gap growth will degradate electrical performance it must be 
controlled and the best way to achieve this control is by 
decreasing the ID/OD ratio. 

3. The innermost ring must be made of iron and is now segmented instead 
of being continuous. This change has been made because the iron ring 
would always reach the limiting strength long before the filamentary 
Wound composite ring(s) reached their strength limit. To overcome 
this limitation a "segmented" inner ring is now proposed for use on 
the magnetically suspended flywheel system. The important point to 
note is that the inner iron ring will have all the necessary magnetic 
properties but will consist of a number of pie shaped segments which 
are bonded to the inside diameter of the first filamentary wound ring. 


The iron ring thus has no stiffness in a "hoop” or tangential 
direction and presents an "inner loading" on the filamentary wound 
composite ring to which it is attached. 

The result of making the above changes in the original GSFC design are all 
intended to bring the original design's excellent theoretical performance to 
fruition in a physical system. 

2. Rotor Analysis Methods 

It is clear that the design dimensions for the GSFC rotor must evolve in 
parallel with the magnetic suspension and motor generator designs (as they 
impact on the dimensions and weight of iron inner ring). Obviously then, the 
most useful rotor design and analysis tools are those which most closely model 
the real physical system and are convenient to apply as the iron inner ring 
design evolves. One such tool has already h$en reported by Kirk and Huntington 
[6, 7, 8, 9] and was used for this project. This tool, an interactive computer 
code called FLYANS (FLYwheel ANalvSis) . represents the most complete stress and 
failure analysis ever performed on a pierced disk multiring rotor. A second 
tool has been developed for this project and is a computer code 

called FLYSlZE (FLYw heel SIZE ) . to size all dimensions for a given rotor 
design. Each of these codes will be discussed in more detail in the following 
sections. 

2.1. FLYANS - Flywheel Analysis Computer Code 
A. Background 

Shown in Appendix 1 is the background theory which was used in developing 
FLYANS. The computer code was originally written by R.A. Huntington [9], as 

part of his Master's thesis, while working with his thesis advisor, Dr. J.A. 

» 

Kirk. The code has since been modified by Dr. Kirk, and the modified version 
was used for this project. 


Shown in Figure n.3 is s schematic diagram of the multiring rotor which is 

modeled in the FLYANS code. Each of the rings are the aame axial thickness and 

♦ 

the stresses in each ring consist of a hoop or tangential atress (o.) and a 

* # tf 

radial stress (o r ). If power is being put in or taken out of the system there 
is an additional shear atress (t^) in each ring. This stress is assumed to be 
zero for this design (i.e., power transfer occurs ’’slowly*') although the FLYANS 
code is set up to allow any stress input the user wishes. It is also assumed 
that the flywheel rings are in a state of plane stress, meaning that there ia 
no variation of the and stresaes in the axial direction. 

The materials which comprise the multiring rotor are modeled as 
homogeneous linearly elastic, specially orthotropic [10] materials, with 
material properties specified in the radial and tangential direction. For this 
project ring #1 is always segmented iron and ring #2 through ring #n are CELION 
6000/EPOXY. The properties for these 2 materials are shown in Table n.l. For 
comparison purposes a number of other materials and their properties are also 
shown. The source of information for properties on Celion 6000/EPOXY was a 
Celanese data sheet provided to us by Mr. Alan Hannibal (Section Manager 
Materials & Process, R & D Department), of Lord Corporation. The other 
materials listed in the table were from an earlier report to GSFC by Dr. J. A. 
Kirk [5]. These other materials are shown here because they constitute the 
working database used by the FLYANS code. It should also be pointed out that 
any new or hypothetical materials can easily be added to the database with 
minimal effort. 

The total stress distribution in one ring of the multiring flywheel is the 
superposition of the five stress distributions due to the following: 



TABLE 3L1 - LIST OF MATERIAL PROPERTIES 
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1. Rotation of the ring at constant angular velocity. 

2. Interaction with adjoining rings due to rotational expansion. 

3. Interference assembly of the rings. 

4. Residual stresses due to curing. 

5. Angular acceleration of the entire assembly. 

The stress dirtribution for the entire flywheel is the summation of the above 5 
stresses for each flywheel ring. 

Of the 5 stress distributions given above, no. 3, interference assembly, 
is under the direct control of the designer. It can therefore be stated that 
it is the goal of the designer to choose an interference pressure set which 
will maximize the flywheel stored energy per unit weight (termed SED). 

B. Interference Stresses 

It will be instructive at this point to consider the hypothetical example 
of how interference stresses interact with rotational stresses in a simple 2 
ring "pierced disk" rotor. 

Shown in Figure n.4 is the stress distribution which occurs when 2 rings 
of the same material are interference assembled (no rotation). When the 
interference stress distribution is added to (i.e., superposed on) the 
rotational stress distribution the net result is as shown in Figure n. 5. In 
Figure n»5 the stresses have been made nondimensional by the factor: 

p 1 «*b* (units are psi) 

where ; 

p x * mass density for the first ring of the assembly (value is 

3 .2 

weight density in Ib/in divided by g * 386 in/sec ) 

w ■ rotational speed (rad/sec) 

b ■ outer radius of the flywhe«l( inches) 
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The solid line shown in each of the plots in Figure n.5 represents the stress 
distribution which occurs when the 2 rings are spun without any interference 

assembly present. The dotted lines show the stress distribution when the 2 

♦ •* 

rings are interference assembled and then spun. Has interference assembly 
actually helped? The answer is yes, because the largest value of 
nondimens ional stress limits performance, and this value has been reduced when 
interference is applied between the rings. Consider the lower plot in Figure 
H. 5. If the working tangential stress of the material, o Q , is constant, then 
the limiting value of b* is ■ 0.97, with no interference present, and 

0.94 with interference present. For a fixed value of b and o Q it is clear the 
interference assembled flywheel has a larger m. Since flywheels are 
characterized by the amount of kinetic energy per unit weight they ‘achieve, the 
interference assembled flywheel will have a larger SED compared to a 
noninterference assembled one. Why? Because the flywheel iB always limited by 
its rotational stresses and £®r a given the interference assembled flywheel 
can rotate faster, and thus have a larger kinetic energy for the same flywheel 
weight. 

C. Constraints 

There are two principal constraints which limit the performance of the GSFC 
multiring flywheel. These are: 

1. The flywheel must not fail. 

2. The flywheel must not exceed a user specified displacement at the 
inner radius (to control the rotor/stator air gap). 

The failure constraint is discussed in detail in Appendix 2. FLYANS allows 
the user to choose between the following 3 failure criterion: 


1. modified TSAI-HILL 

2. maximum STRESS 

3. maximum STRAIN 

and FLYANS vill then return the limiting value of .p x « a b* for the user selected 
geometry and materials. Extensive use of the FLYANS program by Dr. Kirk has 
shown that the TSAI-HILL failure criterion is generally the most conservative 
and it was therefore used for all work on this project. It should be noted 

i 

that a prediction of flywheel performance based on one of the above failure 
theories is a prediction for the maximum or "burst" performance of the rotor. 

The second constraint is on air gap growth and the FLYANS code applies 
this in addition to the failure criterion. Either a failure criterion or 
"inner radius displacement" always limits multiring rotor performance.. 

What about interference? Interference does not limit performance; it 
enhances it. FLYANS analyzes a user specified design assuming no interference 
and then asks if the user wishes to apply interference. If the user elects to 
apply interference* he is asked to either input the ring to ring interference 
pressure he desires or to permit the program to use its own algorithm to choose 
an interference pressure set. The interference pressure selection algorithm is 
documented in Appendix 1 and details may be found there. If the program 
selects the interference pressure set then the user must specify the maximum 
allowable radial interference (i.e., Z INTERFERENCE ■ radial mismatch between 2 
rings divided by the nominal radius of the rings). Typically 0.2 to 0.5Z 
interference is used. The FLYANS code will then compute the "best" pressure 
set and will return the new SED and p x w a b a values. Typically the SED will 
increase 20Z to 200Z over the identical design without any interference. 
Finally* the FLYANS code returns the partial assembly pressures for each of the 


ring*. Rings are always added from the inside out, so knowing the partial 
assembly pressures and ring thickness, it is straightforward to calculate the 
assembly forces (this is exactly what FLYSIZE does). 

D. Summary 

The FLYANS code will analyze multiring flywheels including the presence of 
an iron inner ring. The required input parameters are summarized in Table D.2 
and typical output parameters are shown in Table D.3. The code will select an 
interference pressure set which maximizes the specific energy density (SED) of 
the multiring rotor. When this code was applied to two 1600 Wh multiring 
flywheel designs for this project the SED went from 25 Wh/lb to 41 Wh/lb in one 
case, and from 19 to 42 Wh/lb in a second case. The detailed output for these 
two runs is presented in Section 2.3. 


2. 2. FLYSIZE - Flywheel Sizing Computer Code 

A second computer code has been developed this project 
and the program listing is shown in Appendix 3. The purpose of this code is 
to take output results of FLYANS and combine them with the specifics of a user 
design, in order to produce detailed information which can then be used to 
generate engineering drawings for each of the flywheel rings. The user 
supplied inputs are: 

1. The number of rings in the wheel* 

2. The outside radius of the wheel, (inches) 

3. The inside radius ratios for each of the rings* 

4c The calculated or specified % interference for each interface* 

5. The user specified maximum inner displacement ratio*, (nondimensional, 
equal to actual displacement divided by flywheel outer radius) 

6. The required stored energy, (Wh) 






6 . 


Table II. 2 


FLYANS INPUT PARAMETERS 

1. Number of rings in the flywheel, (dimensionless) 

2. Inner radius ratios of each ring, (dimensionless) 

3. Flywheel materials, each specified by 

a. Name and identification number 

b. Orthotropic ratio, /E 0 /E r , where 

Eg ■ tangential direction modulus of elasticity, (psi) 

E r ■ radial direction modulus of elasticity, (psi) 

c. Tangential modulus of elasticity, (psi) 

d. Weight density, (lbs/in ) 

e. Working stress iu tension in the tangential direction, (psi) 

f. Working stress in compression in the tangential direction, (psi) 

g. Working stress in tension in the radial direction, (psi) 

h. Working stress in compression in the radial direction, (psi) 

i. Shear working stress, (psi) 

j. Material cost, ($/lb) 

k. Residual stress parameter, (dimensionless) 

4. What material each ring is made from (chosen from the user supplied 
database in item 3). 

5. The maximum inner radius displacement ratio (i.e. gap growth divided by 
flywheel outer radius). 

6. The maximum permissible ring interference %• 


I ft b ls II. 3 


FLYAHS OUTPDT 

1. Specific energy density (SED), (Wh/lb). 

2. Volumetric energy density (VED), (Wh/ft 3 ). 

2 2 

3. Limiting value of b , (psi). 

4. Which ring limits performance. 

5. Interference pressures at each ring interface 

a. after complete assembly. 

b. to assemble. 


7. The limiting value or p x «*b*, (pci)* 

8. The calculated SED, (Wh/lb)* 

»• » 

3 

9* The calculated VED, volumetric energy density, (Wh/ft )* 

10. The partial assembly pressures, (ksi)* 

3 

11. The weight densities for each of the flywheel rings*, (lb/in ) 

The FLYSXZE code will then compute the following design values. 

1. The inside and outside radius for each flywheel ring, (inches) 

2. The radial mismatch at each ring interface, (inches) 

3. The flywheel thickness, (inches) 

4. The weight of each ring, (lbs) 

* 

5. The miuimum conical taper which must be present on each ring to allow 
interference assembly, (degrees) 

6. The maximum flywheel rotational speed, (rpm) 

7. The total flywheel weight, (lbs) 

8. The maximum inner radius displacement (air gap growth), (inches) 

9. The partial assembly forces (assuming a 1-1/2° ring taper and 0.1 
interface coefficient of friction), (lbs) 

The calculation formulas which are used to generate the above outputs are 
explained in detail in Appendix 4. Shown in Table n.4 is a summary of the 
design parameters which are utilized or calculated by the FLYSIZE code. The 
symbols shown in this figure are consistent with those used in the FLYSIZE 
program (listed in Appendix 3) and will permit an easier understanding of 
FLYSIZE logic. 


* Taken directly from the FLYANS run. 


f l a(n) 
RR(n-l) 
INF(n) 
;PP(n) 
F(ri) 


SYMBOL 


rU), r(2) — r(n) 

WT(1) , WT(2) WT(n) 

a(l), a(2) a(n+l) 

RR(1), RR(2) RR(n) 

INF(l) , INF(2)— -INF 
PP(1), PP(2)— PP(n) 
F(l) , F(2) F(n) 


SYMBOL TABLE 

INPUT/OUPPUT 


MEANING 


I 

Identification 

I 

Weight density 

0 

Weight 

0 

Radius 

I 

Inside radius ratio 

I : 

X Interference 

I 

Partial assy, pressure 

0 

Partial assy, forces 

I 

Outer radius 

I 

Inner ring disp. 



The FLYSlZE code he* been applied to tvo specific design* for this project 

end the re'sults will be presented in section 2.3. 

» * 

♦ f 


3. Result* 

Both the. FLYANS end FLYSIZE computer codes were used to design two rotors 
for this project, and the design drawings for each' are presented in Appendix 7, 
These two designs are identified as NASA-12 and NASA-20, and they represent the 
results which were obtained from running the FLYANS code over 40 times. 

Shown in Table II. 5 through II. 8 are the FLYANS output for the NASA-12 and 
NASA-20 designs. The results shown in these 4 tables are summarized below. 

1. Table II. 5 shows the results for a 6 ring rotor with ID/OD « 0.5. 

This design is labeled NASA-12. The design is unoptimized and has no 
interference between rings. The SED (specific energy density) is 
25.0 Wh/lb and this number includes the presence of an iron inner 

• ring. Displacement of the inner radius (which is equal to the air gap 
growth due to rotation) is less than .006 times the outer radius of 
the flywheel. 

2. Table II. 6 shows the results for the NASA-12 design with interference 
between rings. The SED has been improved to 41.0 Wh/lb by the effect 
of prestressing (through interference). The inner radius displacement 
is still less than .006 times the outer radius of the flywheel. 

3. Table 11*7 shows the results for a 6 ring rotor with ID/OD «* 0.4. 
this design is labeled NASA-20. The design i6 unoptimized and has no 
interference between rings. The SED is 18.5 Wh/lb and the inner 
radius displacement is less than .004 times the outer radius of the 
flywheel. 

4. Table II. 8 shows the results for the NASA-20 design with interference 
between rings. The SED has been improved to 41.9 Wh/lb and the inner 
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radius displacement is still less than .004 times the outer radius of 
the flywheel. 

One of thie additional, and very informative, outputs from the FLYANS code 
are plots of residual and rotational stresses vs. 'radius ratio for each design 
being considered. These plots are shown in Appendix 6 for the following 
stresses: 

1. Radial residual stress vs. radius ratio 

2. Tangential residual stress vs. radius ratio 

3. Radial, stress (at max. speed) vs. radius ratio 

4. Tangential stress (at max. speed) vs. radius ratio 

and have been obtained for both the unoptimized and optimized NASA-12 and 
NASA-20 designs. 

In order to apply the FLYANS results to a specific design xt is necessary 
to specify the following parameters. 

1. Total stored energy in Wh, a value of 1600 Wh is being used for this 
project. 

2. The outer radius of the flywheel in inches, a value of 10 is being 
used in this project. 

Once these values are specified, the FLYSIZE code can be run using the FLYANS 
output for each design being considered. When FLYSIZE is applied to the 
NASA-12 and NASA-20 designs the results are as shown in Tables II. 9 through 
11.12. Based upon the results shown in these tables it is very convenient for 
the designer to choose the dimensions of each ring of the flywheel. The method 
used to select the ring diameters and tolerances is included in Appendix 5, and 
the detailed drawings for NASA-12 optimized and NASA-20 optimized are shown in 
Appendix 7 of this report. 


NASA-12 UNOPTIMIZED 6 RING CELION 6000/EP0XV FLYWHEEL 


RING 

NO. 

<--) 

IN. R. 
RATIO 
t— ) 

INSIDE 

RADIUS 

(IN) 

RUN ON 12- 
OUTSIDE 
RADIUS 
(IN) 

26-19B3 AT 11 :06s 55 
RADIAL 7. 

MISMATCH I NT 

(IN) tlO A -3) 

RING 

WEIGHT 

(LBS) 

MIN. 

TAPER 

(DEG) 

1 

0.500 

5. 000 

5.200 ‘ 



B. 1 


2 

0.520 

5.200 

6.000 

0.000 

2- 3 0.0 2- Z 

6.9 

0.0 

3 

0.600 

6.000 

7 , 000 

0.000 

3- 4 0.0 3— 4 

10.0 

0.0 

4 

0.700 

7.000 

Q. 000 

0.000 

4- 5 0.0 4- 5 

11.5 

0.0 

5 

0 . BOO 

8.000 

9.000 

0. 000 

5- 6 0.0 5- 6 

1 *5 • 0 

0 . 0 

6 

0.900 

9.000 

10.000 



14.6 



RING THICKNESS (INCHES) . 


4.4 




RPM (ROT/HIN) 



30.B80 




SED (WH/LB) . . 



25.0 




VED(WH/FT A 3) 



1,984.4 




TOTAL STORED 

ENERGY <WH) 


1 . 600 




TOTAL FLYWHEEL WEIGHT (LBS) ... 

64.0 , 




MAX 

. INNER RAD. DISP. ( INCHES) . 

0.060 


* 


PARTIAL ASSEMBLY PRESSURES FORCES 


RING 

3 

ADDED 

TO 

PREVIOUS 

RINGS. .... 

0.0 

kPSI 

0.0 

10*3 

LBS 

RING 

4 

ADDED 

TO 

PREVIOUS 

RINGS 

0.0 

kPSI .... 

0.0 

10" 3 

LBS 

RING 

5 

ADDED 

TO 

PREVIOUS 

RINGS 

0.0 

kPSI .... 

0, 0 

1 0 A 3 

LBS 

RING 

6 

ADDED 

TO 

PREVIOUS 

RINGS 

0.0 

kPSI .... 

o.o 

10'%? 

LBS 


RING WEIGHT DENSITY RING WEIGHT DENSITY 

( — ) (LB/IN A 3) ( — > (LD/IN"3> 


1 0.286 6 0.055 

2 0.055 

3 0. 055 

4 0.055 

5 0.055 


DATA FOR THIS RUN IS IN FILE B: RUN- 12. DAT 


Table II. 9 


ORIGINAL PAGE IS 
OF POOR QUALITY 


NASA-12 OPTIMIZED 6 RING CELION'EPOXY FLYWHEEL, 
RUN ON 12-26-1903 AT U)07r40 


RING 

IN. R. 

INSIDE 

OUTSIDE 

RADIAL 

% 


RING 

MIN, 

NO. 

RATIO 

RADIUS 

RADIUS 

MISMATCH 

I NT 


WEIGHT 

TAPER 

(--) 

< — > 

< IN) 

(IN) 

(IN) 

(10 

5) 

(LBS) 

<DEG> 

l 

0.500 

5.000 

5.200 





4,9 


2 

0.520 

5.200 

6.016 

0.016 

2- 3 3. 0 

r>m» 

3 

4.2 

0 • V 

pt 

v- 

0.600 

6.000 

7.01B 

0.010 

3- 4 3.0 

n- 

4 

6. 1 

0.4 

4 

0.700 

7 . 000 

8.021 

0.021 

4- 5 3.0 

4- 

5 

7.0 

0.4 

55 

0 . 000 

8,000 

9.024 

0.024 

5- 6 3.0 

5- 

6 

7,9 

0.5 

6 

0.900 

9,000 

10.000 





0.9 



RING 

THICKNESS { INCHES) , 


2.7 






RPM (ROT/MIN) , 



39, 531 






SED (WH/LB) . . 



41.0 






VED (WH/FT'N3) , 



CJ 

n 

5 

r« 

n 






TOTAL STORED ENERGY <WH> ...... . 1,600 

TOTAL FLYWHEEL WEIGHT (LBS) . . . 3-9.0 

MAX. INNER RAD. DISP. (INCHES) . 0.060 


PARTIAL ASSEMBLY PRESSURES FORCES 


RING 

•Mf 

ADDED 

TO 

PREVIOUS 

RINGS. , 

• • • 

3.9 

fsPSl 

43.7 

10-3 

LBS 

RING 

4 

ADDED 

3 

PREVIOUS 

RINGS. . 

• < • 

4.6 

IPS I .... 

58.7 

10-3 

LBS 

RING 

5 

ADDED 

r 0 

PREVIOUS 

RINGS. . 

• * * 

4.5 

KPS I 

66.7 

10-3 

LBS 

RING 

& 

ADDED 

TD 

PREVIOUS 

RINGS. . 

• • * 

4.2 

KPSI .... 

72. 1 

10-3 

LBS 


RING WEIGHT DEN5ITY RING WEIGHT DENSITY 

(--> (LB/IN a 3) 4 — « (LB/ IN' -3) 


1 0. 286 6 0.053 

2 0.055 

3 0.055 

4 0.055 

5 0.055 

DATA FOR THIS RUN IS IN FILE 3jRUN-12~0.DAT 


Table 11.10 


ORIGINAL PAGE 19 
OF POOR QUALITY 


RING 

NO. 


J 

2 

3 

4 

* 5 

6 

original page is 

OF POOR QUAWTY 


NASA-20 UNOPTOMIZED 6 RING CELION/EPOXY 6 RING FLYWHEEL 




RUN ON 12' 

-24-19B3 AT 

09» 53i 06 



IN. R. 

INSIDE 

OUTSIDE 

RADIAL 

7. 

RING 

MIN. 

RATIO 

RADIUS 

RADIUS 

MISMATCH 

INT 

WEIGHT 

TAPER 

(— > 

(IN) 

(IN) 

(IN) 

<10"~3) 

(LBS) 

(DEG) 

0.400 

4.000 

4.200 



8.1 


0.420 

4.200 

5.000 

0.000 2- 

3 0.0 2- 3 

7.0 

0.0 

0.500 

5.000 

6.500 

0.000 3- 

4 0.0 3- 4 

16.4 

0.0 

0.650 

6.500 

8.000 

0.000 4- 

5 0.0 4- 5 

20.7 

0.0 

0. 800 

8.000 

9.000 

0.000 5- 

6 0.0 5- 6 

16.2 

0.0 

0.900 

9.000 

10.000 



18. 1 



RING THICKNESS (INCHES) 3.5 

RPM (ROT/M IN) .................. 27, 4B4 

SED (WH/LB) 10.5 

VED (WH/FT'*3) 1,596.7 

TOTAL STORED ENERGY (WH) ....... 1,600 

TOTAL FLYWHEEL WEIGHT (LBS) ... B6.6 

MAX. INNER RAD. DISP. < INCHES) . 0,040 


PARTIAL ASSEMBLY PRESSURES FORCES 


RING 

3 

ADDED 

TO PREVIOUS 

RINGS. 

... . 0.0 

kPST . . . , 

, 0.0 

lO'-S 

LBS 

RING 

4 

ADDED 

TO PREVIOUS 

RINGS, 

0 * * • 0# u 

kPSI . . . , 

, 0.0 

10 A 3 

LBS 

RING 

5 

ADDED 

TO PREVIOUS 

RINGS. 

. . . . 0,0 

kPSI 

■ 0.0 

10*3 

LBS 

RING 

6 

ADDED 

TO PREVIOUS 

RINGS. 

.... 0.0 

kPSI .... 

, 0.0 

L0"3 

LBS 

RING 

WEIGHT 

DENSITY 


RING 

WEIGHT 

DENSITY 



( — > 


(LB/ 

IN A 3> 


C— > 

<LB/IN A 3) 




1 

r> 


4 

5 


0.286 
0 . 055 
0.055 
0.055 
0.055 



DATA FOR THIS RUN IS IN FILE B: RUN-20. DAT 


Table II. 11 


NASA-20 OPTIMIZED 6 RING CELIQN/EPOXY FLYWHEEL 


RUN ON 12-26-1983 AT 09i02il6 


RING 

IN. R. 

INSIDE 

OUTSIDE 

RADIAL 7. 

RING 

MIN. 

NO. 

RATIO 

RADIUS 

RADIUS 

MISMATCH INT 

WEIGHT 

TAPER 

<— > 

<--> 

(IN) 

(IN) 

(IN) 

( 10“~3) 

(LBS) 

(DEG) 

1 

0. 400 

4 , 000 

4.200 



3.6 


2 

0.420 

4.200 

5.025 

0.025 

2-3 5.0 2- 3 

3. 1 

O. 6 

3 

0.500 

5.000 

6. 533 

0.033 

3-4 5.0 3- 4 

7.2 

o.a 

4 

0. 650 

6.500 

B. 040 

0.040 

4- 5 5.0 4- 5 

9.1 

0.9 

5 

0.800 

8,000 

9.045 

0.045 

5-6 5.0 5- 6 

7. 1 

1.1 

6 

0. 900 

9.000 

10.000 



B. 0 



RING 

THICKNESS t INCHES) . . 


2,4 




RPM(ROTZMIN) , 


* * V * « • 

41,372 




SED (WH/LB) . . , 


• » » ** • 

41.9 




VED (WH/FT 3) , 



3,613.1 




TOTAL STORED 

ENERGY <WH> . 

1 • • » i * 

i , 600 




TOTAL FLYWHEEL WEIGHT (LBS) ... 

30.2 




MAX. 

INNER RAD. DISP. (INCHES) . 

0.040 

1 





PARTIAL ASSEMBLY 

PRESSURES FORCES 


RING 

3 

ADDED TO PREVIOUS 

RINGS. . , 

B.5 

kPSI .... 01.5 

lO'-S 

LBS 

RING 

4 

ADDED TO PREVIOUS 

RINGS. . , 

9.9 
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4. Discussion 

i * 

For a given fo*® factor (XD/OD) rotor there are 2 principal advantagea of 
uaing multiring construction versus winding the rotor aa one continuoua ring* 

Theae are: 

* 

1. The residual atreaaea due to fabrication are known to be amaller for 
•mailer size ringa. 

2. The individual ringa of a nultiring design can be interference 
aaaenbled. The condition will significantly improve performance and is 
much better than building in rotor pres tress through controlled winding 
tension. 

The main disadvantage of multiring construction is that each ring interference 
must be bonded together upon assembly. This binding must be sufficiently 
strong so that the interference radial stresses (due to rotation) will not cause 
the outer rings to separate from the inner ones. Since the required strength of 
the interface bond is always less than (or equal to) the transverse strength of 
the graphite/epoxy, the same epoxy used in ring fabrication will accomplish the 
task of bonding the rings together. It can be concluded that the advantages of 
multiring fabrication outweigh the disadvantages and that the multiring 
construction with interference assembly represents the best manufacturing 
method or the NASA-12 optimized and NASA-20 optimized designs. . 

The NASA- 12 optimized design will produce 41.0 Wh/lb at its maximum 
rotational speed of 39,581 rpm. This design has an inside diameter of 10 
inches an outside diameter of 20 inches, and an axial thickness of 2.7 inches. 
It will store 1600 Wh of energy when it is running at maximum speed, and its 
air gap growth will not exceed .060 inches at 39,581 rpm. The design has an 
iron inner ring weighing 4.9 lbs and the total weight of the flywheel is 39.0 


lb«. There are 6 ring* in this design and rings % through 6 are Ceiion 
6000/EPOXY (filamentry wound composite). There is light interference (.3%) 
between ail Celion 6000/EPOXY rings and the assembly forces to press the rings 
together do not exceed 75,000 lbs. 

The NASA-20 optimised design will produce 41.9 Wh/lb at its maximum 

* 

rotational speed of 41,372 rpm. This design has an inside diameter of 8 
inches, an outside diameter of 20 inches, and an axial thickness of 2.4 inches. 
It will store 1600 Wh of energy when it is running at its maximum speed, and 
its air gap growth will not exceed .040 inches at 41,372 rpm. The design has 
an iron inner ring weighing 3.6 lbs and the total weight of the flywheel is 
38.2 lbs. There are 6 rings in this design and rings 2 through 6 are Celion 
6000/EPOXY. There is heavy interference (.5%) between all Celion 6000/EPOXY 
rings but the assembly forces to press the rings together do uot exceed 125,000 
lbs. 

It was generally observed that to control growth ot the air gap it is 
desirable to decrease the inside diameter of the rotor. To control the 
stresses which result from this action it is necessary to increase ring 
prestressing through an increase in % interference between rings. To 
illustrate this point, had the NASA-20 design been limited to 0.3% interference 
between rings, it would develop only 35.2 Wh/lb versus the 41.9 Wh/lb which it 
will develop with 0.5% interference between rings. 

It must be emphasized that the design values presented for NASA-12 and 
NASA-20 are for the flywheel pushed to its limiting failure speed. If the 
flywheel is operated at .707 of this limiting speed the maximum flywheel 
•tresses will be 1/2 their failure values and the SED will be 1/2 of its 
limiting value. If it is required that the usable stored energy be 1600 Wh and 
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the flywheel ie assumed to cycle between .707 of its limiting speed end .25 of 
its limiting speed, then the calculated flywheel thickness must be increased, 
in both the, NASA-12 and NASA-20 designs, by a factor of 2.3. Obviously, then, 
r>en a flywheel is said to store 1600 Wh, it could only deliver this amount of 
energy ifi 

1. The flywheel is cycled from its maximum rotational speed to zero speed 

2. The flywheel did not fail at its maximum rotational speed. 

Clearly, the final rotor design for a magnetically suspended flywheel must 
proceed with consideration for the speed range of the motor/generator and the 
iron ring weight requirements for the magnetic suspension and motor/generator 
magnetic paths. The FLYANS/FLYSIZK tools used in this project are the best way 
to revise the rotor design as other system components evolve. 
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III. FABRICATION/ASSEMBLY 

A. Materials. 

Two candidate graphite epoxy materials exist*for this application. These 
materials are; 

giber 

Celanese Celion 6f> n 0* or 17,000 
Union Carbide Thornel 3C0, WYP6 1/0 

The fiber should be without twist and have a matrix compatible sizing. 
Epoxy 

Shell EPON 826* with Jeffamine Chemical Company D230 Curative 
NARMCO 5213 

* Recommended for this task. 

B. Techniques 
1. Fabrication 

Residual transverse and/or circumferential stresses should be held to a 
minimum. Prestressing using fiber tension has not proven practical. Either 
precorapression of thin rings (via interference assembly) or a flexible urethane 
resin appears at present to be the best approach. Urethane resin use would 
require further study and testing before a practical matrix could be specified. 

A controlled environment, such as room temperature and a relative humidity 
less than 50%, are desirable during winding. These conditions control 
variations in material properties. In addition, to assure quality materials, 
the resin should be stored carefully according to the manufacturer's 
recommendations and shelf life strictly adhered to. To verify the resin, a pad 
should be cast with each flywheel for subsequent testing. 


Initial and final tie-off procedures as recommended in Lawrence Livermore 
Laboratory report UCRL-15437 should be adhered to wherever appropriate. 

The cure and post cure cycles should be designed to optimize the 
properties of the resin and to minimize residual'stresses. Usually, a heat-up 
cycle at a controlled rate to cure temperature followed by a cool down cycle 
also at a controlled rate can achieve both of these goals. Report UCRL-15437 
page 11 thoroughly discusses this. 

Fiber volume fraction (nominally 60-65%) should be controlled to within a 
±2X tolerance. The void content should be maintained at less than 3%. These 
quantities should be checked at the inside, middle and the outside of each 
ring. The surface finish should be free of bubbles, pock marks, resin rich 
areas or exposed fibers. 

2. Assembly 

All composite rings should be tested both ultrasonically and 
radiographically before assembly. Only high quality rings should be included 
in a multiring configuration. Subsequent flaw detection should also be 
administered in order to assure no fracturing occurred during assembly. 

Assembly procedures as specified on the assembly drawing must be strictly 
followed. Each ring pair should be static balanced for mating assembly and 
overall static balance will be obtained by removing material from the inner 
ring only. 

C. Drawings 

The drawing package is found in Appendix 7 and contains drawings for two 
different ID/OD ratios. The sponsor will include the appropriate package for 
the design selected to be built. 


D. Candidate Fabricators 

The demise of the Department of Energy program in Mechanical Energy 
Storage Technology (MEST) has reduced the number of possible candidate flywheel 
fabricators to only two. They are: •• 

Lord Corporation (A.J. Hannibal) 

1635 West 12tb Street 
P. 0. Box 10039 
Erie, PA 16514-0039 

and 

Hercules Corporation (P. Ward Hill) 

Aerospace Division 
P.O. Box 98 
Magna, Utah 84044 

At present, the strongest candidate from the standpoints of technical 
capability, interest and corporate commitment is the Lord Corporation. 

E. Testing 

Over the past several years two experimental facilities have been active 
in testing composite material flywheels. These are: 

1. The Oak Ridge Flywheel Evaluation Laboratory, located at the Oak Ridge 
Y-12 plant near Knoxville, Tennessee. 

2. The Applied Physics Laboratory of Johns Hopkins University in Laurel, 
Maryland. 

Both facilities have been funded by the Department of Energy and both have been 
significantly dismantled as a result of DOE no longer funding a flywheel 
effort. In addition, both facilities were built with the sole intention of 
testing shaft driven flywheels. Therefore, a major portion of their 


experimental equipment consisted of a high speed air turbine and air supply 
source, neither of which is required for the GSFC design. 

Our curr.ent assessment of GSFC spin testing requirements would indicate 
that the following experimental equipment would be' required to evaluate a GSFC 
magnetically suspended flywheel design: 

1. An air tight enclosure having a volume of at least 7 ft 3 (3 feet in 
diameter by 1 foot high). 

2. A containment ring to surround the flywheel system. 

3. One or more optical viewing ports to permit cameras to view the 
enclosure interior. 

4. Instrumentation to measure pressure, temperature, wheel runout, wheel 

t 

wobble, etc. 

5. An air flushing system to exhaust epoxy and composite material fumeB 
from the enclosure. 

6. A capture device which will restrain the spinning rotor so it does not 
destroy the suspension or motor/generator electronics if the rotor 
decouples. 

_2 

7. Vacuum equipment to pump the enclosure to 10 TORR or better. 

8. Remote viewing and data collection capabilities. 

9. Controllers for the motor/generator electronics to either run at a 
constant speed or cycle between two reference speeds for fatigue 
testing. 

Our best estimate is that while some of the facilities of the two laboratories 
may still be available and applicable to GSFC requirements, it is not cost 
effective to adapt them for GSFC testing. 


IV. SPECIFICATION FOR COMPOSITE FLYWHEEL RING-ROTOR 
Introduction . 

The objective of this task is to fabricate, assemble and static test a 
multiring rotor for the NASA/GSFC "Spacecraft Flywheel Power System" program. 
Previous studies by NASA/GSFC showed the merits of considering Inertial Energy 
Storage as a possible viable alternative to electrochemical or fuel cell power 
storage devices for long life, high cycle rate, extended use, space 
applications. These studies showed that a magnetically suspended ring rotor, 
made of filament wound construction, containing an integral ironless armature 
motor-generator rotor* operating in an evacuated hermetically sealed enclosure 
should be an efficient means for space energy storage use. This effort 
addresses the first systems element, the "Filament Wound" rotor. 

Scope 

The scope of this phase of the program includes fabrication, assembly and 
static testing of the filament wound multiring rotor. Analysis has established 
the best fiber/epoxy materials and the method of prestressing of multiring 
elements to maximize energy storage capacity. 

Requirements 
A. Size 
1. Assembly 

OD - 20.000" +0.010" 

ID- 8.000" +0.010" 
or 10.000" +0.010" 

W - 2.700" +0.005" 

or 2.400" +0.005" 


2. Components 

Rings No. l f 1A 

See drawings No. 7-1 & 7-1A in Appendix 7 
Rings No. If 2A *' 

See drawings No. 7-2 & 7-2A " " 

Rings No. 3, 3A 

See drawings No. 7-3 & 7-3A 11 " 

Rings No. 4, 4A 

See drawings No. 7-4 & 7-4A " " 

Rings No. 5, 5A 

See drawings No. 7-5 & 7-5A " v 

Rings No. 6, 6a 

See drawings No. 7-6 & 7-6A ” " 

f 

Note: 

Mandrel cylindricity control is critical. Each mandrel will have the 
taper stated on the appropriate drawing with a 152 tolerance allowable. This 
must be strictly adhered to in order to meet the assembly load requirements. 

B. Material Selection 
Fiber 

Celanese Celion 6000 graphite fiber will be used. 

Eposcy 

Shell EPON 826 with Jeffamine Chemical Company D230 curative will be used. 
Mixing will be as prescribed by the manufacturer. 


C. Processing 

Wyap Lav -Up 

1 . Tensioning 

Fiber tensioning should be controlled to minimize residual transverse 
and/or circumferential stresses. It is desirable to strive for a zero residual 
stress winding technique. 

2. T emperature 

Winding will take place in a clean "dust-free" area with room temperature 
controlled to 75°F ,+10’F. Relative humidity will be controlled to be less than 
50*. 

3. Resin 

Resin will be stored carefully according to the manufacturer's 
recommendations and shelf life strictly adhered to. To verify resin, a pad 
should be cast with each mix for subsequent testing. 

4. Fiber Tie-Off 

Initial and final tie-off procedures as recommended in UCRL-15437 will be 
adhered to. 

D. Cure Method 

The cure method and post cure cycle shall be designed to optimize the 
properties of the resin and to minimize residual stresses. A heat-up cycle at 
a controlled rate bo 230°F followed by a cool down cycle also at a controlled 
rate can achieve both of these goals. See page 11 of UCRL-15437 for control 
rates. An oven with a good commercial temperature controller will be adequate 
for this operation. 


E. Finishing 

The fiber volume fraction shall be controlled to within a tolerance of 
±22. Void content shall be maintained at less than 3%. The surface finish 
shall be free of bubbles, pock marks, resin rich areas or exposed fibers. The 
machined outside diameter of each ring shall be smooth and free of loose fiber 
ends. 

F. Assembly 

Drawing Nos, 7-7 and 7-7A (Appendix 7) cover assembly of the ring-rotor. 
Assembly procedure stated must be followed implicitly. Only composite rings 
that passed both ultrasonic and radiographic inspection will be used. Ring 
alignment must be controlled during the "nesting" process. Press anvil and ram 
faces must be parallel and smooth to prevent surface damage of the rings. 

After each ring is assembled, the sub-assembly will be tested ultrasonically 
and radiographically and certified before moving on to the next assembly 
procedure. The completed assembly will be tested and certified (see Section 
4). 

G. Testing 

1. Fiber 

A certificate of physical and mechanical properties shall be included with 
each spool of graphite fiber. Only materials that exceed the minimum 
manufacterer 's requirements will be used. 

2. Ep ox y 

Pads taken from each batch mixed (see Section III - Resin) will be tested 
in accordance with the manufacturer's recommendations to insure that 
conformance to specifications has been met. 


3. Rings 

Fiber volume fraction and void content aha.ll be checked at the inside, 
middle t>nd outside of each ring. 

4. Ring Assembly 

All composite rings shall be tested both ultrasonically and 
radiographically before assembly. Only high quality rings shall be included 
for multiring assembly. Subsequent flaw detection chall be administered during 
the subassembly build-up to assure that no fracturing occurred during the 
assembly process. 

H. Balancing 

Prior to assembly of rings each ring should be static balanced and 
oriented so, upon assembly, any imblaance compensates that which exists in 
other rings. After complete assembly final static balance will be achieved by 
removal of material from the inside of the inner ring. Final static balance 


will be 300 micro grams or less 


V. SUMMARY 


The design definition of a composite fiber ring rotor has been completed. 
Tvo specific designs have been included with this ‘report, both capable of 
storing a maximum energy of 1600 Wh. Both designs account for the presence of 
a iron inner ring and both utilize Celion 6000/EP0XY as the composite material 
for the fiber composite rings. 

The first design is identified as NASA-12 and consists of a total of 6 
rings with the flywheel having an ID/OD ratio of 0.5. The first ring is iron 
and the remaining 5 rings are interference assembled Celion 6000/EPOXY. The 
projected maximum Specific Energy Density (SED), for the interference assembled 
flywheel, is 41 Wh/lb (90.2 Wh/kg) and 25 Wh/lb (55.0 Wh/kg) for the identical 
design assembled without interference. The air gap growth, from zero to full 
•peed, will not exceed .060 inches. 

The second design is identified as NASA-20 and also consists of a total of 
6 rings, with the flywheel having an ID/OD ratio of 0.4. The first ring is 
iron and the remaining 5 rings are interference assembled Celion 6000/EPOXY. 

The projected maximum SED is 42 Wh/lb (92.4 Wh/kg), for the interference 
assembled flywheel, and 19 Wh/lb (41.8 Wh/kg) for the identical, geometry 
assembled without interference. The air gap growth, from zero to full speed, 
will not exceed .040 inches. 

The major differences between NASA-12 and NASA-20 are in the level of 
interference which is built into the ring. NASA-12 is designed for 0.3% 
maximum ring to ring interference while NASA-20 has 0.5%. NASA- 20 has a 
smaller ID/OD ratio and thus a smaller maximum air gap growth. Both designs 
appear equally viable but before either is built it is recommended that two of 


the Celion 6000/EPOXY rings, from the KASA-20 design, be fabricsted, strain 
gaged, and pressed together. Details of this recommendation may be found in 
Section VI. ■ ' 

In addition to producing 2 specific designs, this project has also 
demonstrated the usefulness of two computer codes, FLYANS (FLYvheel ANalySis) 
and FLYSIZE (FLYvheel SIZE), for efficiently carrying out multiring rotor 
design ^ including the presence of an iron inner ring. Since the rotor design 
may have to be modified as the motor/generator and magnetic suspension evolves, 
the computer codes will allow this to be done as easily as possible. In 
addition, both codes are available for use in projecting the types of future 
composite material properties which will most benefit the GSFC magnetically 
suspended flywheel. 

A review of the spin test requirements for the GSFC eastern indicates 
existing exerimental facilities will not be applicable. Two of the most active 
facilities, Oak Ridge and Johns Hopkins University Applied Physics Laboratory, 
have been effectively dismanteled due to their loss of DOE funding. Even if 
these facilities were fully functional, the unique capabilities of the GSFC 
system (i.e. to spin its own rotor) effectively eliminate the need for a major 
item in these facilities, namely the high speed air turbine. It would thus 
appear that a new spin test facility, which takes advantage of the exerience 
gained from the DOE program, would be most cost effective for NASA/GSFC. 


VI, RECOMMENDATIONS 

r. 

The results presented in this report show that interference assembly of 

. *» 

rotor rings will be required to achieve competitive energy densities with shaft 
driven flywheel designs. It is suggested that 2 rings (Number 5 and Number 6) 
of the NASA~20 design be built, instrumented with strain gages, and then 
pressed together. The experimental results can be compared to the theoretical 
predictions in order to enhance our understanding of pressing assembly 
performance. After these 2 rings have been assembled, they can be 
disassembled, and ring No. 5 machined to generate ,3Z interference (making it 
identical to ring Number 5 in the NASA-12 design). The rings can then be 
instrumented and pressed together again, and results can be compared to 
theoretical predictions. Finally, the rings can be disassembled and then 
radially cut in order to determine the residual stresses present in each ring. 
This test sequence requires only 2 rings and will provide the great** 3 t 
information return for the coat invested. 

In characterizing the composite ring-rotor, the iron ring was treated as a 
uniform profile mass load. Before a practical rotating assembly can be built, 
the iron ring must be defined in detail to identify the magnetic suspension 
portion of the assembly. In addition, the inner composite ring must be further 
defined to accommodate the ironless armature portion of the motor/generator. 
This also includes provisions for the magnet assembly. 

Both of these refinements will necessitate recalculating the stress 
patterns and analytically determining the true operating characteristics of 
this design. While it is felt that the new calculations will not seriously 
de-rate the ring-rotor^ defined in this study, one must know what the actual 
©gsrating capacity of this design is before test hardware is built. 


Completion of this follow-on effort will result in a detailed definition 
of the rotating assembly. This approach is a necessary and logical sequence in 
pursuit of a total energy storage system design. 
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APPENDIX 2 


FAILURE ANALYSIS 


from 


Huntington, R. A., M. S. Thesis, "Stress 
Analysis and Maximization of Performance 
for a Multiring Flywheel," University of 
Maryland, Mechanical Engineering 
Department, Jan. 1978. Thesis advisor: 


To determine the performance of the multiring flywheel, it has been stated 
that a -suitable failure criterion must be applied to the stress distribution. 
This appendix will present several failure criteria, including the maximum 
stress, maximum strain, Tsai-Hill, and Tsai-Wu failure theories [1, 2-5]*. 

Their development, applications, and limitations will also be presented along 
with their specific application to the multiring flywheel. 

A simple test for the determination of a failure stress is a uniaxial 
tensile test in which a specimen is continuously loaded until rupture. The 
results of this test may be applied directly to a design situation in which 
uniaxial laods along a known direction are found. However, virtually all real 
situations involve stress states which are biaxial or triaxial with both normal 
and shear stress components. Since it would be very impractical, if not 
impossible, to perform a failure test for every combination of stress which may 
be foreseen, it becomes necessary to develop a failure criterion, based on a 
limited number of experiments, to predict if any arbitrary stress State will 
cause material failure. This prediction of failure stress has been studied in 
isotropic materials since the ninteenth century. Two widely used criteria, the 
Treska and von Mises, have emerged from these studies. Both assume that 
failure is induced solely by the action of shear stresses. The von Mises 
criterion further defines the distortional energy - the strain energy not 
associated with dilatation - as a parameter which reaches a critical value at 
failure. Modifications of both of these criteria (for isotropic materials) 
have been proposed to eliminate original restrictions [6-8]. Also, a 


* Denotes references at back of this Appendix. 


generalization of the von Mises criterion has been proposed by Hill [3] for 
slightly anisotropic metals. However, the materials considered for the 
multiring flywheel consist of both isotropic metals and orthotropic metals and 
orthotropic composites. Therefore, either one very general or several 
individual failure criteria are necessary to meet the requirements of the 
multiring flywheel. Two of the failure theories to be discussed, the Tsai-Hill 
and Tsai-Wu criteria, meet this requirement of generality since they both 
reduce to the von Mises failure theory for isotropic materials. The other 
failure criteria to be discussed should not be applied to isotropic materials, 
although it may be applied to pseudo-isotropic composite laminates. 

The complexities of failure prediction for a composite material is easily 
portrayed by a brief discussion of the stress states in an isotropic versus an 
orthotropic materials under an arbitrary loading. For an isotropic material, 
the state of stress may be transformed at will to any convenient reference 
plane (usually the principle stress or maximum shear directions) without change 
in material properties. However, if the stress state of the composite is 
transformed, both the elastic properties and material strengths vary. While 
the variations of elastic properties follow the transformtion of a fourth rank 
tensor, the transformation of strength is not well known. Therefore, a failure 
criterion for a composite material must relate all components of stress to a 
plane of symmetry in which strength properties are well known. This adds to 
the number of parameters necessary to characterize a failure surface. A 
failure theory for composites is further complicated by the variety of modes of 
failure which exist. The three major modes are fiber fracture, matrix cracking 
and delamination [9]. Therefore, it is necessary that the criterion predict 


•tress limits for several failure modes. However, it is not necessary to 
actually predict the mode of failure which will occur under a given loading. 

Maximum Stress Criterion 

The maximum stress failure criterion assumes that failure will occur in a 
composite material when any one of the stress components exceeds a critical 
value 1-2]. Each component of stress is considered independently from the 
others and all of the critical stress values are independent. This failure 
criterion is given by the following equation set: 

( 2 , 1 ) 

are the uniaxial compressive, 

If any of the inequalities given by Eqn. 2.1 are violated, failure is 
postulated to occur. This is a very simple criteria to use and mainly because 
of this, it is often used to predict failure of composite materials. There 
are, however, several major restrictions on the general application of this 
method. The first of these requires that the components of stress be 
transformed to the material symmetry axes. But for the multiring flywheel, the 
stress distributions are given in the r-0 plane while this is also the assumed 
plane of material symmetry. Therefore, ths restriction does not limit 
application for the multiring. The most important restriction on the 
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application of this criterion is simply is accuracy. This is demonstrated by 
Figure 2.1, which is a plot of predicted and experimentally determined failure 
stresses versus the angle between the line of action of a uniaxial load and the 
fiber orientation of the composite material. It should be noted that the 
loading axis uses a logarithmic scale so that the apparently small deviations 

i 

represent large errors in the criterion. It should also be noted that this 
failure criterion predicts discontinuities in the failure strength, associated 
with changing failure modes, while the experimental data shows smooth 
variations. 


Maximum Strain Criterion 

The maximum strain failure criterion is very similar to the maximum stress 
criterion in application. The only major difference between these criteria is 
the assumption that a critical value of strain limits material strength rather 
than a critical value of stress (as with the maximum stress criterion). These 
strain values may be transformed to stress terms by application of Hooke's Law 
yielding the following equation set which describes this failure criterion in 
the case of plane stress: 


‘ 'r < ° r " °e < °r 
- o= < o 9 - v re o r < °e 


r6 r 

“ T r0 < T r0 < T r0 


( 2 . 2 ) 


where: etc. have been previously defined. 


Comparison of Eqns. 2.1 and 2.2 6hows that the only difference between 
these failure criteria is the poisson effects between the normal stress 
components. Along with their similarities, these criteria share identical 
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Figure 2.1 
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off-axis strength (after ref. 2 ) 


restrictions (i.e.: stress components must be transformed to material symmetry 

axes and inaccurate predictions of failure stress). The latter of these is 
ahovn by' Figure 2.2 which is a plot of predicted and experimentally determined 

i » 

failure stress similar to Figure 2.1. Again large errors are seen in the 
predicted failure stress over a wide range of loading orientation. 

Tsai-Wu Criterion 

The preceding failure criteria each have several inadequacies which impair 
their agreement with experimental data and their application. An obvious 
method of improving data fit is to increase the number of independent 
parameters in the criterion. The next failure criterion to be discussed, the 
Tsai-Wu theory [5, 10, 11], achieves this with a tensorial failure theory which 
relates linear and quadratic stress components. The basic postulate is that 
there exists a failure surface with an equation of the form: 

F i a i + F ij a i a j “ 1 (2.3) 

where: F^ * 2nd rank tensor of strength parameters 
F^. •* 4th rank tensor of strength parameters 

These strength tensors will obey known tensor transformations for 
coordinate changes and are asusmed to be symmetric. The former demonstrates 
that it is now possible to determine the strength properties of an anisotropic 
material in an arbitrary plane once the terms of the strength tensor are 
established. This is a major advantage of this failure theory although this 
feature is not needed to analyze multiring flywheel failure. It may be noted 


that the linear stress terms allow for variations in the tensile and 
compreasive strengths while the quadratic terms define a cloaed ellipaoidal 
failure surface which intersects each stress axis provided that the following 
stability relations Isl are satisfied 


F ij ± F ii F jj 

where: repeated indices are NOT summations 


(2.4) 


The failure theory, as represented thus far, is applicable to general 
anisotropic materials with triaxial loading ; however, the multiring flywheel 
deals with specially orthotropic materials under plane stress loadings. For 
this case, Eqn. 2.3 may be expanded, yielding: 


F l B r + F 2 a B * F 6 T r0 


+ F Ur 


+ *12 Ve + F 22°e 


|'I» p T* 

*66 T re 


(2.5) 


Using Eqn. 2.5, the F's may now be evaluated with experimental data. All 
but Fj 2 may be found by uniaxial tensile and compressive testing and assuming 
that the positive and negative pure shear strengths are equal, the F's become: 


*11“ 

T c 

V * 

r r 

a a 
r r 

F 1 

a T o c 
r r 


1 


_ c _ T 

a e a e 

P 22 

T r* 

F 2“ 

°e °e 

Q 

> 

CD 
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( r^> 2 

T re 
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F 66 ~ 

F 6 
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The stress interaction term, must be determined with a combined 

stress loading arrangement. A great deal of work has been carried out to 
determine optimal tests for this parameter [5, 10, 11}; however, conclusive 


results have not been presented in these references. This identifies a 
significant disadvantage of this failure theory since information vhich may 
determine F^ is not generally available for the composite materials being 
considered for the multiring flywheel. To examine the effect of variations of 
on the predicted performance of a flywheel rotor, is defined using 
Eqns. 2.4 and 2.6 as: 

f 12 " n ® ( 2 « 7 ) 

V °r a r °6 a 0 

where: -1 < C < 1 for stability conditions 

Figure 2.3 shows the effect of variations of F^» within these stability 
limits, on the predicted SED of a single ring flywheel (a^ ■ .75). Since the 
maximum prediction is three times greater than the minimum, it is obvious that 
the SED, and other performance parameters, is very sensitive to change in F^« 

Despite this disadvantage, the Tsai-Wu failure, the Tsai-Wu failure theory 
has good agreement with experimental data if F^ is properly selected. This is 
demonstrated by Figure 2.4 which shows predicted and experimental data for 
off-axis strength tests on a boron-epoxy composite. Despite this improvement 
in accuracy over the previously discussed failure theories, the lack of 
information needed to determine prevents the further use of this failure. 
However, an additional failure criterion, the Tsai-Hill theory [1, 2—4) 9 which 
is similar to the Tsai-Wu theory, will be discussed next. 

Tsai-Hill Criterion 

The Tsai-Hill criterion is a generalization of the von Mises yield 
criterion for isotropic materials. It was first proposed by Hill in 1948 [3] 



o 



for slightly anisotropic metals and latter adapted C 196t») to composite 
materials by Tsai (4] • As described in these references, the criterion 
contains only quadratic stress terms. However, linear terms have been added to 
allow for variations in the tensile and compressive strengths. The resulting 
criterion, which will be termed the modified Tsai-Hill criterion, is identical 


to Eqn. 2.5 with the F's defined by Eqn. 2.6. However, this criterion 
contrasts the Tsai-Wu in that the interaction term is not independent of 
the other strength parameters. For the materials considered in this project, 
®12 - - T ir Therefore, the parameter C from Eqn. 2.7 becomes: 


C 


1 

1 

\ 


a T a c 
r r 

„ T_ c 

a e °6 


( 2 . 8 ) 


For filamentary wound composite or isotropic materials, it is obvious that 
C is within the stability requirements given by Eqn. 2.7. The agreement of the 
modified Tsai-Hill failure criterion with experimental data is shown in Figure 
2. A. Also, the value of C given by Eqn. 2.8 for the previously cited single 
ring flywheel has been indicated on Figure 2.3. It is easily seen that the 
agreement with experiment is not as good as is seen for the Tsai-Wu theory. 
However, this sacrifice of some accuracy is necessary to eliminate the problem 
of determining Fj^ for the Tsai-Wu theory. Therefore, the modified Tsai-Hill 
criterion will be used instead of the Tsai-Wu theory for the multiring 
flywheel. However, the development of the remaining theory will be general 
enough to allow a simple transition back to the Tsai-Wu criterion if F^ values 
become available in the future for the materials considered in this study. 
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FLYSIZE LISTING 


(This program if written for an IBM PC using BASIC.) 


ORIGINAL PAGE IS 
OF POOR QUALITY. 


Hi’ a FILYS I ZE 


3B i~u S 
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Page 1 


10 


-FLYS 1 ZE 




20 DIM A (15), F(10) , RMM (15) , GAMMA (15) , LABEL* < 15) 

30 D 3: M PROMPTS* (5) , PROMPT? * ( 1 0 ) , RR ( 1 S j , I NF (IS) , PP (15) v GAMMA* (10), PH I D I < 1 0 ) 
40 SCREEN 2 : KEY ‘OFF 
50 WIDTH 80 
60 CLS 

70 •' INPUT USER DATA— 

80 60 SUB 5000 


90 J 
100 
110 
120 
130 
140 
1 50 
160 
3. 70 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 


? «.**#** 




NOW DO COMPUTATIONS" 

CALCULATE RING THICKNESS- 
PI "3 » 14159 

T * 1 728 # KE / ( P I % B'”2 # OLD ) 

■’ — CALCULATE WHEEL RPM — 

RPM» 1 07 . 6 3, # ( BCR ( RHOOMEGABSG / GAMMA ( 1 ) ) ) / B 

? - - CALCULATE- MOM I IMAL R I IMG RAD I US VALUES- 

FOR I i*l TO N-l-1 
A (I) « RR(I)#B 
NEXT I 

CALCULATE PARTIAL ASSEMBLY FORCES- 
AND MINIMUM TAPER ANGLE FDR EACH RING PAIR- 

MUb. 1 
PH ID“1 - 5 
CON V-P 1 / 1 80 ! 

PHIR»PHID*GONV 
FOR I «2 TO N— 1 

F ( I ) «2#P I *A(IJ mPF < I ) t ( MU+TAN ( PH I R ) ) 

PH ID I ( I ) « ATM ( INF ( I ) *. 001 JKA ( I > /T> /COIMV 




280 NEXT I 

290 •’ 

300 FOR 
310 RMM ( I ) 


JALCULATE RADIAL M I S MATCH- 
TO N-l 

-■ INF (I) 1' A (I) # .001 


ii ihi im i»« w» win j»i«.fim i«ii i.n . 


•30 RMM (0) -•••0 
140 R'MM ( N ) «0 

•50 -CALCULATE RING 


OUGHTS- 


L ; 


360 

370 

380 

390 

400 

410 

420 

430 

44 0 

450 

460 

470 

480 

490 

500 

510 ! 

5000 

5005 

5010 


THE WEIGHT DENSITIES FOR THE RINGS ARE GAMMA ( 1 ) THRU GAMMA (N) 

’THE WEIGHT DENSITIES ARE READ IN AS USER DATA 






.III, H'H III* lllll |HI< »»*>» lit. 


WTOTAI. 
FOR I 
RING AREA 
RING VOL 
WT ( I ) 

W TOTAL - 
NEXT i: 


- 0 


. TO N 

« PI * (A(M-l) -"-2 » A (I) ■•"•2) 
* T * RTNGAREA 
« RING VOL # GAM'-A ( I ) 

W TOTAL WT ( I ) 


-END OF COMPUTATIONAL LOOP 


* MU) »>•(♦ 


REM 

REM 

REM 

’ OUTPUT THE 

GO SUB 7000 
END 


sBULTS- 


•INPUT & FILE CREATION SUBROUTINE 

•INPUT IS FROM THE KEYBOARD OR FROM A USER SPEC IF IE FILE- — - 

•THE USER FILE IS B ii RUN-X X X X ,, DAT AND MUST BE ALREADY AVAILABLE ON 


Os F 7 L_ V S X Z; El - IB tf”‘t S 


ORIGINAL PAGE IB 
OP POOR QUALITY 

pr i nt ed on 1 2-26- 1 983 at U * 1 9 s 30 


15015 
5020 
5025 
5030 
5035 
5040 
5045 
5050 
5055 
5060 
5065 
,5070 
! 5075 
! 5080 
, 5085 
I 5090 
1 5095 
51.00 
5105 
5110 
'5115 
1 5120 

! ki 1 o nr 
| wJ 4 . *iW w 

*5130 
,5135 
15140 
5145 
5150 
5155 
5160 
5165 
.5170 
1 51 75 
‘51 80 
5185 
15190 
15195 
5200 
jS205 
1 52 10 
5215 
,5220 

M A w 

5230 

5235 

15240 

,5245 

5250 


>'*a *m* i** »»» •in 


-DRIVE HO IF THE INPUT 18 SPECIFIED FROM THE KEYBOARD A USER 
-FILE OF THE SAME FORM WILL BE CREATED ON DRIVE EL AND THE INPUT 
-DATA WILL BE WRITTEN THERE » 

CLS . 

LOCATE 2,5 

PRINT "WILL YOU INPUT FROM <F>ILE" 

LOCATE 3,5 

PRINT ” OR CIOEYBOARD" 

LOCATE 5, 5 
PRINT "ENTER F OR I 
A#~ I NKE Y $ 

IF A#«" " GOTO 5065 




ii . 


IF A$«»K" OR Af>- 


f ii u u 


THEN FLAG 7- 1 : GOTO 5090 


IF A#-"F" OR A#'»"f" THEN FLAG7~2s GOTO 5090 
GOTO 3065 

IF FLAG 7.-2 GOTO 5105 
LOCATE 7,5 

PRINT "A FILE WILL BE CREATED FOR YOUR INPUT DATA ON DRIVE Bs " 

LOCATE 8,5 

PRINT "THE NAME OF FILE IS Bs RUN-XXXX „ DAT" 

LOCATE 9,5 

PRINT "XXXX SHOULD BE AN INTEGER BETWEEN 0 AND 9999" 

LOCATE U,5 

INPUT;; "ENTER XXXX NOW " ID* 

F I LENAME#- " B s RUN- " • * - ID# + " , DAT " 

IF OR THEN OPEN FILENAME# FOR OUTPUT AS #2 

IF A$«»F" OR A$*"f" THEN OPEN FILENAME# FOR INPUT AS #1 
PROMPT'D# * "ENTER A RUN TITLE UP TO 80 CHARACTERS LONG" 

PROMPT 1# * "NUMBER OF RINGS" 

PROMPTS# « "ENTER RADIUS RATIO FOR EACH RING INTERFACE" 

PROMPTS# <0> * "INSIDE RADIUS RATIO" 

PROMPTS# < 1 5 » " OUTER RAD I US RAT I 0 C MUST- 1 - 0 ) " 

PR0MPT4# * "ENTER INTERFERENCE FOR EACH RING INTERFACE" 

PROMPTS# * "ENTER THE PARTIAL ASSEMBLY PRESSURES" 

PROMPT?# < 1 ) * "FLYWHEEL OUTER RADIUS ( W> " 

PROMPT?# (2) » "STORED ENERGY < WH) " 

PROMP 1 7# C 3 ) 1 == " RHO 1 < OMEGA #S > # #2 ( PS I ) " 

PROMPT?# < 4 ) * " Si ED ( WH / LB ) " 

PROMPT?# (5) « "VED<WH/FT##3> " 

PROMPT?# < 6 ) » " IMOIM-D I M . D I SP . L I M I T < 1 0"""-3 > " 

’ — - — -INPUT AN BO CHARACTER (MAX ) DESCRIPTIVE RUN TITLE— 

LOCATE 14,5 * PRINT PRQMPTO# 

LOCATE '15,1 i IF FLAG7-1 THEN INPUT RUNT I TLE# ELSE INPUT #1, RUNT ITLE# 
IF FLAG7-1 THEM WRITE #2, RUNT ITLE# 

f»i f-i 


*• ttminn h.i 




a 4 


60 


i: 

WilwWv 


,, U/ »l.. 


70 


INPUT NUMBER OF RINGS- 
LOCATE 1,36 s PRINT PROMPT 1# 

LOCATE 1,51 ,s IF FLAG7-1 THEN INPUT N 
IF FLAG7-- 1 THEN WRITE #2,N 

? - -INPUT RING RADIUS RATIOS- 

LOCATE 2,10 SPRINT PROMPT 2# 

FOR I * l TO N 


..BE INPUT #1, N 


ORIGINAL PAGE ® 

i of poor quality 
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'5283 

5290 

3295 

5300 

5305 

5310 

6313 

5320 

5325 

5330 

5335 

3340 

5345 

5350 

5355 

5360 

5365 

6370 

5375 

5380 


IF m THEN PRINT USING "RING ### TO #tttt "s 1-1:1 
LOCATE 1*2,51 ' ’ 

IF FLAG7.k1 THEN INPUT RRM) ELSE INPUT # 1 H RR ( I ) 

IF FLAG 7. *1 THEN WRITE #2,RR<I> 

NEXT I’ 

INPUT OUTER RADIUS RATIO, MUST EQUAL 1.0- 

LOCATE I +2 , 20 : PR I NT PROMPTS* < 1 ) 

LOCATE 1+2,51 

IF FLAG7=1 THEM INPUT RR (N+l ) ELSE INPUT #l„RR(N+3> 

IF FLAG7=;I, THEM WRITE #2,RR(N+1) 

INPUT RING TO RING INTERFERENCE 

LOCATE 14,8 : PRINT PR0MPT4* 

FOR I ~2 TO N-l 
LOCATE 1+14,1 

rWljrt» , NB '[ INTERFERENCE RING ### TO RING #### (10'"-3 IN/ IN) “s I ; 1 + 1 
LOCATE 1+14,51 ' 

IF FLAG7=1 THEM INPUT I NF ( I ) ELSE INPUT #1,INF<D 
IF FLAG7=1 THEM WRITE #2, INF (I > 

NEXT I 
m q 


530E 
5390 
5395 
5400 
5405 
5410 
5415 
5420 
5425 
5430 
5435 
3440 
5445 
5450 
5455 
. 5460 
5465 
5470 
5475 
5480 
5485 
5490 
5495 
5500 
5305 
5310 
ITY E 
5515 
5520 

w v«a.|L \„f 

5530 

nsi&'TKr, 

WWWW 

5540 

5545 


* INPUT MI SC),, VALUES* 

FOR I = i TO 6 

LOCATE 1+4,25 : PRINT PR0MFT7* (I ) 

s rP FUA(3% ” 3 ‘ THEW * NPUT VAR a) ELSE INPUT #1, VAR CD 
IP I- LAB7= 1 THEN WRITE #2, VAR (I) 

NEXT I 
B = VAR < 1 ) 

KE * VAR (2) 

RHOOMEBABSQ * VAR (3) 

SED * VAR (4) 

VET.) = VAR (5) 

D I SP » VAR ( 6 ) *BY9 , 999999E-04 

CLS 

" INPUT PARTIAL ASSEMBLY PRESSURES- — — 

LOCATE 1,14 : PRINT PROMPTS* 

FOR I =2 TO M~ 1 
LOCATE 1+1,20 

." PP R:I:NS ## # ADDED TO ALL PRIOR RINGS (KSI ) " j; I + :i, 

LuCA 1 b , { 63 

IF FLAG7.k3. THEN INPUT PP(I) ELSE INPUT #1„ PR (I) 

IF FLAG7=:L THEM WRITE #2, PPM) 

NEXT I 

INPUT THE RING WEIGHT DENSITIES IN LB/1W'73 ...... 

CLS 

L. 0 L A 1 b o , 2 0 s I"'' R j. N I 1 ' b |\| T E R T H E W E I G H T D E l\t S I T Y 0 F E A C H R I N (3 1 1 
WFP T R 1,4,3 I|PRINT " IF PRESENT RtWG & ALL FOLLOWING RINGS HAVE THE SAME 

FOR 1=1 TO N 
LOCATE 1+5,5 

F49INT USING "ENTER DENSITY FOR RING ### (L..B/IM A 3) " n 1 1 
IF FLAG 7=1 THEN INPUT GAMMA* (I ) ELSE INPUT #1 , GAMMA $ < I ) 

IF F LAG 7= 1 THEN WRITE- #2, GAMMA* <1 > 

1 STAR » I 

IF GAMMA* ( I ) = "S" OR GAMMA* ( I ) ="s" THEN 5565 


DENS 
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VAL(GAMMA*(I> ) 


Page 4 
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5550 GAMMA (I) = 

15555 NEXT I. 

5560 GOTO 5585 

5565 ? SETS -ALL SUBSEQUENT RINGS TO SAME DENSITY IF S ENTERED- 

5570 FOR I “I STAR TO N 

5575 GAMMA ( I ) =VAL < GAMMA* ( I STAR- 1 > ) 

5580 NEXT I 

5585 LOCATE 1+6,10 ;■ PRINT " * >K>!< ^ ^ H< Jtc>K ^ >K >KEND OF USER INPUTKW******)}!** " 
5590 CLOSE #1 

5595 IF FLAG’X— 1 THEN CLOSE #2 
5600 RETURN 

7000 OUTPUT SUBROUTINE 

7005 LABEL* ( 1 ) « 11 R I NG TH I CKNESS ( I INCHES ) " 

70 1 0 LABEL* ( 2 ) « » RPM ( RQT/M I. IN) " 

70 1 5 LABEL* ( 3 ) - " BED < WH/LB) " 

7020 LABEL* (4) -"VED < WH/FT--3) " 

7025 LABEL* ( 5 > * " TOTAL STORED ENERGY ( WH ) " 

7030 LABEL* (6) -•"TOTAL FLYWHEEL WEIGHT (LBS) " 

: 7035 LABEL* < 7 ) * " MAX . I INNER RAD . D I SP . ( I INCHES ) " 

. 7040 CLS 
7045 LOCATE 8,14 

r 1 7050 PRINT "ENTER THE FIRST LETTER FOR ONE OF THE FOLLOWING CHOICES" 

| 7055 LOCATE 12,23 

1 7060 PRINT "<S>CREEM OUTPUT TO SCREEN" 

.7065 LOCATE 14,23 

7070 PRINT " < E >IND » „ „ „ * . END THE PROGRAM " 

7075 LOCATE 17,27 

7080 PRINT "YOUR CHOICE" ; 

'7085 A* ~ INKEY* 

7090 IF A*™ "8" OR A*k»"b" GOTO 7105 
'7095 IF A*“"E" OR A*“"e" THEN END 
’7100 IF A**" " GOTO 7085 

7.1.05 OPEN "SCR IN I! " FOR OUTPUT AS #1 : GOTO 71.10 
'7110 CLS 

,71 15 LOCATE 3. , 40 -LEM ( RUNT I TLE* ) /2 
7.1.20 PRINT RUNT I TLE* 


,7125 

LOCATE 2,27 

7,1,30 

PRINT #1, "RUIN 

■* y | •-.jMJij 

LOCATE 3,3 

7140 

PRINT #1 , "RING 

"RADIAL" TAB (56) "X. 

-7:1.45 

LOCATE 4,4 

7150 

PRINT #1, "NO,, " 


- 7 140 PR I NT # 1 , " R I NG " TAB ( 9 ) " I IN „ R , " TAB ( 1 El ) " I INS I DE " TAB ( 27 ) " OUTS I DE » TAB ( 40 ) 

’ " RAD I AL " TAB ( 56 ) " X. " TAB ( 64 ) » R I NG " TAB ( 75 ) " M I IN " 

7145 LOCATE 4,4 

7150 PRINT #1), "NO., " TAB (10) "RATIO" TAB (IS) "RADIUS" TAB (27) "RADIUS" TAB (39) » 
* 1 SNATCH " TAB ( 55 ) " I NT " TAB ( 63 ) " WE I GHT" TAB ( 74 ) " TAPER 11 
t 7 155 LOCATE 5,3 

'7160 PRINT #1. !, " (— ) " TAB ( 10) " ( ~) " TAB ( 19) " ( IN) " TAB (28) " ( IN) " TAB (41 ) " ( IN) 

:TAB ( S3 ) " ( 1 0 3 ) " TAB ( 64 ) " (LBS ) " TAB ( 74 ) " ( DEG ) " 


7165 

LOCATE 

6, 1 


'7170 

PRINT 1 

H , STRING* (80,45) 

M 

JJ 

17175 

K ”= 6 



"71.80 

FOR I > 

* 1 TO IN 


...7185 

LOCATE 

I +!< , 3 : PRINT #1, 

USING "###"? I 

7190 

LOCATE 

I-H<’l0 SPRINT #1 , 

USING "#,, ###" ; RR ( I ) 

3 7 195 

LOCATE 

I +l< <i 18 SPRINT #1, 

USING "###.###": A(I) 
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V LC)C A 1 E 3. +K , 27 : PR I NT # 1 , US X NG '*'### „ ### " : A < I + 1 ) +RMM ( I ) 

5 IF I«1 THEN GOTO 7230 

0 LOCATE I +-K , 39 : X F KN THEN PRINT #1, USING "#.###•': RHN { X ) 

5 LOCATE. I+K r 45s IF I<N THEN PRINT # 1 , USING "# -##" I 1*1 

0 LOCATE I+K, 52: IF KN THEN PRINT #1, USING "#7# ,I ;INF < I > 

5 LOCATE I+Kj, S6: IP* KN THEN PRINT #1, USING""# -##" s I „ H-l 

-> LOCATE I +K , 75; IF I < N THEN PR I NT # 1 , US I NG " #7# " ; PH I D I ( I ) 

^ LOCATE I+K, 63 ' 

PRINT #1, USING "####* #" ; WT ( I ) 

:> NEXT I 
5 L * 6 

:> LOCATE J.*L, 14: K*1 : GOSIJB 10000 

5 LOCATE K-L, 45: PRINT #1, USING "###„#" T ORIC 

1 LOCATE I *L-I- 1 p 1 4 : !<“2 : GOSUB 10000 ’ OF 1 

3 LOCATE K-L+lp 45s PRINT #1 * USING "######, . " : RPM 

) LOCATE I H-L+Sji 1 4 s K*3 : GOSUB 10000 

3 LOCATE I+L+2, 45: PRINT #1, USING "*##.#" ; SED 

] LOCATE K-L-K3, 14: K-4: GOSUB 10000 

3 LOCATE I -i-L+3 » 45 : PR I NT # 1 , US I NG "##### „ , # " - VED 

! I-OCATE I+L+4, 14 s K*5s GOSUB 10000 

5 LOCATE H-L+4,, 45: PRINT # 3. > ( USING "######, „ 11 jKE 

> LOCATE 14 : K«6s GOSUB 10000 

1 LOCATE I*L+5,45 » PRINT #1, USING "####. : WTOTAI 

i LOCATE K-L *6, 14 :K-7: GOSUB 10000 

J LOCATE 1 H-L+ 6 . 1 45 s PR I NT # 1 US I MG " tt „ ### » $ D I SP 

i GOSUB 11000 

i CLS 

' LOCATE 1,15 

i PRINT "PARTIAL ASSEMBLY PRESSURES 

' M~0 

FOR 1*2 TO N-l 
LOCATE K-M, 4 

CAPTIONS “"RING ### ADDED TO PREVIOUS RINGS" 

PR I NT US I NG CART 1 0N$ I + 1 « 

PR I NT STR I NB$ ( 40-L..EN ( CART 1 0N$ ) , 46 ) ? 

LOCATE I +M ,, 4 1 : PR I NT tJS I IMG " #### . # " ; pp ( j ) 

LOCATE K-M,, 49: PRINT "kPSI" 

LOCATE K-M , S3 : PRINT USING " #####.# 10-\3 LBS" s F (I ) # 

NEXT I ■' 

I STAR-- K-M 
CO LX* 1 :: RCJW%»0 

C APT I ON :1$“" RING WEIGHT DENSITY" 

C A P T 1 0 N 2 " ( — ) ( l p / t i\i •■••• p ) 1 1 

LOCATE I STAR't'l , CO LX 
PRINT CART I ON 1$ 

LOCATE ISTAR+2, COLX 
PRINT CART I 0N2$ 

LOCATE l‘STAR+3 , COL’/, s PR I NT STR I NG* ( 25 , 45 ) s 
IF N >5 THEN LOCATE ISTAR+.1. , 40: PRINT CAPTION 1 $ 

IF N >5 THEN LOCATE ISTAR+2, 40: PRINT CAPTIONS^ 

.1.1- N>5 1 HEN LOCATE ISTAE+3, 40: PRINT STRINGS <25„ 45) 

FOR 1-1 TO N ? V 

IF I >5 THEN COL 7, '- :i 40 s R0W7-5 
LOCATE I STAR+3+ 1 --ROWX C0L7, 


FORCE 


WEIGHT DENSITY" 
(LB /I N"'"3 ) " 


I 


B i a F=*L-YE> I ZEE! » IS^SB 
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7475 PRINT USING ” ###"« I 
7 4 BO LOCATE I STAR+3+ 1 -ROW’/. , COL.7. +15 

74S5 PRINT USING "#.###" ; GAMMA ( I > ORIGINAL PAGE IS 

7490 NEXT I ■ OF POOR QUALITY 

7495 IF N >5 THEN ENDREF*/* 3 IBTAR+ 10 
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APPENDIX A 
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4-1 Calculation of Taper Annie 

Shown in Figure 4.1 iB a schematic diagram showing the cross section of 
two rings , Ring (I) and Ring (1+1), which are to'be assembled using pressing 
force F(l). In order to permit these 2 rings to be press assembled, it is 

i 

necessary that point A lie inside (i.e. towards the center line) of point B . 

To calculate the minimum value of 0 which will permit this, assume that point A 
and point B are directly above each other, and further assume that the 
following information is known for these 2 rings: 

T - ring axial thickness, (inches) 

INF(l) * X interference between ring (I) and (1+1), (dimensionless) 

A(I) - nominal radius of the interface between ring (I) and (1*1) * 

( inches). 

It then follows from the definition of X interference that: 

RMM(l) « INF(I) * A(l) 
where : 

RMM(I) «* radial mismatch between Ring (I) and (1+1). 

Then, it can be seen that 0(1), the taper angle, should be selected so that: 
0(1) > TAN -1 (RMM(I)/T) (4-1) 

FLYSIZE calculates the taper angle for each interfering ring pair using the 
equality in equation (4-1) and provides the results in a printed design report 
for the subject flywheel. In general the angle 0(1) is less than 1.5 degrees 
for up to 0.5Z interference. 
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4-2 Calculation of Assembly Far cat 


Shown in Figure 4.2 ia a schematic diagram of the pressing forces (F(I)) 

» * 

which acts on Ring (I). Since the Rings (I) and (1+1) are tapered the pressing 

force will be a maximum at the point where the two rings are flush with each 

other. Assume that at this point the following information is known: 

P(I) ■ partial assembly pressure between Ring (l) and (i+l), (psi) 

p - 0.1. coefficient of friction between the two rings. This value is 

taken as a constant for FLYSIZE and assumes the rings are put 

together using epoxy as an effective lubricant. 

t(I) - p * P(I), interfacial shear stress between the two rings, (psi). 

Since the forces are in equilibrium it is known that: 

1 Forces ■ 0 
z 

Therefore : 

F(I) - t(D * AREA(I) cos 0 + P(I) AREA(I) sin 0 

- AREA(I) (t(I) cos 0 + P(I) sin 0) (4-2) 

where: 

2 

AREA(I) ■ surface area of the interface between the two rings, (in ) 

0 ■ 1.5°, assumed to be constant for all ringss, (degrees) 

The expression for AREA(I) is obtained assuming that the line segment T - A is 
revolved about the center line of the ring. Therefore: 

AREA(I) - 2 * n * (T/cos 0) * (A(l) + RMM(l)/2), (in 2 ) (4-3) 

In practice RMM(l) is at least 2 orders of magnitude less than A(l), so 
equation (4-3) simplifies to: 

A(I) ■ 2 * n * (T/cos 0) * A(l), (in 2 ) ( 4—4) 

If equation (4-4) is substituted into equation ( 4-2), then: 



Schematic diagram of the pressing force 
and interface stresses acting on Ring (I) 




(4-5) 


F(I) - 2 * n * (T/coa 0) * AU) * (p*P(I) cos 0 + P(I)*sin 0) 

* 

F(I) - 2 * n * A(I) * T * P(I) * (0 + tan 0), (lbs) 

FLYSXZE calculates the assembly force between Rings (I) and (1+1) using 
equation (4-5). The actual pressures (P(I)) used by FLYSIZE are the partial 
assembly pressures required to add Ring (1+1) to all previously assembled 
rings* This means that the force calculated using equation (4-5) will be the 
partial assembly force required to add Ring (1+1) to all the previously 
assembled rings. Note that multiring assembly always starts from the inside 
ring outward, and that the first 2 rings which are press assembled are rings 2 
and 3. Ring 1 is always segmented iron and requires no interference assembly. 
Practically this means that it is added last and simply bonded into place. 
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4-3 Calcuation of Ring Thickness 

To calculate the ring thickness (T), assume that the following information 
is known: 

KE ■ the required total stored energy in the flywheel, (Wh) 

3 

VED "the volumetric energy density of the subject flywheel, (Wh/ft ) 

b ■ flywheel outer radius, (inches) 

Since VED ■ KE/(nb^T) by definition then: 

T - 1728 * KE/(n * b 2 * VED), (inches) (4 -6) 

3 3 

where the constant 1728 converts ft to in . 

FLYSIZE calculates the ring thickness using equation (4H>>. 


4-4 Calculation of Maximum Wheel RPM 


To calculate the maximum wheel rpm assume that the following information 
is known: 

p a «*b* ■ a known value from the FLYANS code, (psi); call this K 

p x - mass density of the ring (l);equal to weight density (Y(D) of 
ring (1) divided by g (386 in/sec ) 
b * outer radius of the flywheel (inches) 

Since w is the angular rotational speed (radians/sec), it is equal to the rpm 
(N) as: 

. - sup <«> u - 7 > 

If equation (4-7) is applied to p x « a b* * K then: 


simplifying 


N - 


(ML- rnr~ 

<2)(nXb) J K/Pa 


» - a i fn k) <y W TT) 


„ 187.61 PT~ , , 

M - b y^cTT ’ (rpn) 


(4-8) 


(4-9) 


(4-10) 


FLYSIZE calculates the flywheel rpm using equation (4-10), 


4-5 Calculation of Miscellaneous Values 


To calculate the nominal ring radius (A(I)) values, assume the following 
information is known: 

RR(I) ■ Inside radius ratios for Ring (I), (dimensionless) 
b ■ outside radius of flywheel, (inches) 

Then: 

A(l) » RR(I) * b (4-11) 

FLYSIZE computes the inside radius of each ring using equation (4-11). If 
there is no interference between rings the outside radius of Ring (I) is equal 
to the inside radius of Ring (1+1), namely A(I+1). If there is interference 
then the outside radius of Ring (I) is equal- to the nominal inside radius of 
the next Ring (A(I+1)) plus the radial mismatch (RMM(I)) between Ring (I) and 
(1+1). Recall that RMM(I) has previously been discussed in section 4-1. 

To calculate the weight of each ring assume that the following information 
is known: 

A(I) = nominal inside radius of ring (I), (inches) 

A(I+1) ■ nominal outside radius of ring (I), (inches) 

3 

y(I) « weight density for ring (I), (lb/ in ) 

Then: 

WT(I) = RINGVOL * y(I) (4-12) 

where : 

RINGVOL = volume of ring (I), RING AREA * T, (in 3 ) 

RINGAREA » n * (A(l+1) 2 - A(l) 2 ) 


T ■ ring thickness, (inches) 


FLYSIZE calculates the Weight of each ring using equation (4-12), and then 

adds up all the ring weights to obtain the total weight of the rotor. 

To calculate the maximum inner radius displacement (i.e. gap growth) of 

the flywheel, assume that the following information is known: 

b * f'lywheel outer radius, (inches) 

IRDR ■ inner radius displacement ratio, (dimensionless) 

Since : 

IRDR - *g 
b 

Then : 

gg - (lRDR)*b (4-13) 

where : 

gg * air gap growth or displacement of inner flywheel radius from 0 speed 
to maximum speed, (inches) 

FLYSIZE calculates the maximum air gap growth using equation (4-13). 
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Determining Ring Dimensions Using 
FLYSIZE Results 


As an example of the procedure used to select ring dimensions and 
tolerances * .consider the design NASA-12, rings 3 and 4, as previously shovn in 
Table H. 10. The dimensions of the two rings are summarised below: 


Ring 

Inside 

Outside 

Radial 


Radius 

Radius 

Mismatch 

<-) 

(inches) 

(inches) 

(inches) 

3 

6.000 

7.018 ^ 

0.018 

4 

7.000 

8.021 ^ 



First the radius of each ring is converted to a diameter and a tolerance 
of ±>001 inches is assigned to the diameter. Thus: 

Q.D. (ring 3) » 14.036 ,+ .001 inches 
I.D. (ring 4) ■ 14.000 .+ *001 inches 

Using these specifications the maximum and minimum radial interference at the 
3-4 interface can be obtained and is: 

max. radial interference B 0.019 inches 
min. radial interference = 0.017 inches 
This varia£!m> ip interference in considered acceptable and will not degrade 
perform?*#*!#* It is now necessary to choose a tolerance on the ring taper 
angles so that the design interference will be acceptable after assembly is 
completed. 

Consider the drawing shown in Figure 5.1, showing the interface between 
rings 3 and 4. The nominal taper angle, 0(nom), must be selected so that point 
B on ring 3 lies to the inside of point A' on ring 4. A valve of 0(nom) = 1° 
will accomplish this for all rings in the NASA-12 design so it will be used. 

It is mandatory that the tolerances on the interface taper angle, 0, be 
Selected so that the interference between the top of ring 3 (point A) and the 


14.036 *"[+. 0005 ] 



Figure 5.1 

Ring 3 and Ring 4 interface in the NASA design 




top of ring 4 (point A') be "close to" the interference between the bottom of 
ring 3 (point B) and the bottom of ring 4 (point B'). This means that the 
nominal value of dimension Y must not vary by more than .+.0005 inches between 
ring 3 and ring 4. Assuming the 0(nom) is 1 degree* then: 

Y(nom) "2.7 * tan(l°) ■ 0.0471 inches 

Thus: 

Y(max) - 0.0471 + .0005 - .0476 
Y(min) * 0,0471 - .0005 ■ .0466 

and: 

Tan 0(max) « Y(max)/2.7 
0(max) * 1.01 
Tan 0(min) “ Y(min)/2.7 
0(min) ■ 0.99 

The tolerance on 0 is therefore .+0.01 degrees. 

Practically, this means that the interface between rings 3 and 4 must be 
measured and then machined to assure that the diameter of the top of ring 3 is 
within ,+0.0005 of the diameter of the top of ring 4. And, the diameter of the 
bottom of ring 3 is within j+0.0005 of the bottom of ring 4. 


APPENDIX 6 


Plots of Stress vs. Radius Ratio for 
NASA-12 and NASA-20 Rotor Designs 
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APPENDIX 7 


Engineering end Assembly Drawings for 
NASA-12 and NASA-20 
6 ring flywheel rotor designs 


Drawing 7-1 


OF POOR QOAUTY 


« © £ © *, 

_ 8 2 » t 

? 2 •« « 

•o o Cfl - 

0) CD — c (0 

co JZ ^ to „ 

(0 


os i- w £ r w 

e «>S = o’ 

■-E ge 2 


CO 


O « «' 


x: — co a> q> — 
H d « -d *o £ 
• ■ 

CM 


CM 


V- 

—u - 

o 


o 

¥ 

«! 

< 

E 

-> 

o 


V 


o 

■M 

« 

%m 

• 

c 


• 

o 

a 

• 


c 

• 


< 

o 




r 


r 


CM+| 

zzzza 


O +| 


ms 



Drawing 7-2 
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Drawing 7-3 
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Drawing 7-4 





Drawing 7-5 





Drawing 7-6 



Draftsman 







Note: 1. Coat all Interfaces 
with epoxy prior 
to assembly. . 










Drawing 7 -2 A 





Drawing 7-3A 
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Drawing 7-6A 







Note: 1. Coat all Flag interfaces 
with epoxy prior to 





